Creation of a homogeneous plasma column by means of 
hohlraum radiation for ion- stopping measurements 



o 

< 

in 

(N 



CZ2 

'Hh' 

O 

• »— i 

Or 



Steffen Faik a •*, Anna Tauschwitz ab , Mikhail M. Basko cd , Joachim A. Maruhn ac , Olga Rosmej ce , Tim Rienecker 3 , 

Vladimir G. Novikov d , Alexander S. Grushin d 

"Goethe-Universitdt Frankfurt am Main, Max-von-Laue-Str. 1, 60438 Frankfurt am Main, Germany 

b Helmholtz International Center for FAIR (HICfor FAIR), Max-von-Laue-Str. I, 60438 Frankfurt am Main, Germany 

1 ExtreMe Matter Institute (EMMI), GSI Helmholtzzentrum fur Schwerionenforschung GmbH, Planckstr. I, 64291 Darmstadt, Germany 

d Keldysh Institute of Applied Mathematics (KIAM), Miusskaya sq. 4, 125047 Moscow, Russia 

' GSI Helmholtzzentrum fur Schwerionenforschung GmbH, Planckstr. 1, 64291 Darmstadt, Germany 



Abstract 

In this work, we present the results of two-dimensional radiation-hydrodynamics simulations of a hohlraum target 
whose outgoing radiation is used to produce a homogeneously ionized carbon plasma for ion-beam stopping measure- 
ments. The cylindrical hohlraum with gold walls is heated by a frequency-doubled {A/ = 526.5 fim) 1.4 ns long laser 
pulse with the total energy of E; = 180 J. At the laser spot, the peak matter and radiation temperatures of, respectively, 
T ~ 380 eV and T r as 120 eV are observed. X-rays from the hohlraum heat the attached carbon foam with a mean 
density of pc — 2 mg/cm 3 to a temperature of T ^25 eV. The simulation shows that the carbon ionization degree 
(Z w 3.75) and its column density stay relatively stable (within variations of about ±7%) long enough to conduct 
the ion-stopping measurements. Also, it is found that a special attention should be paid to the shock wave, emerging 
from the x-ray heated copper support plate, which at later times may significantly distort the carbon column density 
traversed by the fast ions. 

Keywords: 2D radiation hydrodynamics, hohlraum radiation and spectra, creation of homogeneous plasma 
conditions, ion stopping in dense plasmas 
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1. Introduction 

Today, the heavy-ion stopping in matter at normal 
conditions yj] is a rather well understood phenomenon. 
The combination of a high-power Petawatt laser facil- 
ity and a large-scale accelerator for heavy ions at GSILj 
offers the unique opportunity to extend this knowledge 
to dense plasmas at high temperatures. Corresponding 
experiments with laser-generated plasmas ]2| |3j |fjj |5|] 
are of crucial importance for the indirect drive scenario 
of heavy ion fusion J6J] and for the ion-driven fast igni- 
tion concept 13,181], but inevitably the biggest challenge 
for those experiments is always to avoid spatial non- 
uniformities in the plasma layer within a sufficiently 
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long lifetime — on the order of a few nanoseconds — 
to measure the stopping power for a bunch of fast ions. 

One straightforward way to form a uniformly ionized 
plasma layer with a constant column density might be 
to heat a planar foil target — as usually being used for 
the measurements with cold matter — by intense direct 
laser radiation. Although at a first glance spatial non- 
uniformities, resulting from the strongly non-uniform 
intensity distribution across the laser focal spot, seem to 
be a formidable obstacle, recent theoretical research [jgl 
confirmed by experiments lUOi ll 111 has shown that un- 
der appropriate conditions sufficiently uniform plasma 
states, suitable for ion-stopping measurements, can nev- 
ertheless be achieved in the direct irradiation scheme. 

Another attractive way to generate a uniform plasma 
state is to heat a sample indirectly with the radiation of 
a millimeter-scale hohlraum. The hohlraum — a cavity 
usually made of a high-Z material, which provides high 
diffusive resistivity for thermal x-rays — thereby is di- 
rectly heated by intense laser pulses to x-ray tempera- 
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tures of tens and hundreds of electronvolts 11121 Il3lll4ll . 
The sample may be either placed inside or close to such 
a hohlraum, usually in the form of a low-density foam. 
Then the inertial and thermal confinement of the uni- 
formly heated high-temperature plasma can be guaran- 
teed either by the hohlraum itself or by an additional 
enclosure for a limited period of time. 

Theoretical modeling of a hohlraum target is a chal- 
lenging task for computational physics since it com- 
bines multidimensional hydrodynamic simulations with 
the solution of the spectral transfer equation for thermal 
radiation. In 2011, 2D (two-dimensional) simulations 
of two hohlraum configurations, one of them already 
used in experiments with the NHELIX and PHELIX 
lasers at GSI II 1511 . were reported II 1611 . In the present 
work, we present the results of new 2D simulations of 
another hohlraum target that is actually used in exper- 
iments with the PHELIX laser and the UNILAC ion 



accelerator at GSI II 1711 . The configuration consists of 
a simple cylindrical hohlraum with gold walls and an 
empty interior. One hole of the hohlraum serves as the 
laser beam entrance. At the other hole, a low-density 
CHO foam inside a copper holder is attached to measure 
the ion stopping in a partially ionized carbon plasma at 
moderate temperatures of T ~ 20 - 30 eV. 

This paper is organized as follows. In Section [2] 
the employed radiation-hydrodynamics code RALEF- 
2D together with the equation of state and spectral opac- 
ities are briefly described. Section [3]gives an overview 
of the target geometry used in the experiments and of 
the corresponding numerical setup for the simulations, 
which have been performed in two steps. In the first 
step, in Section|4j the heating and radiation of the empty 
hohlraum without the foam sample are analyzed. In the 
second step, in Section |5J the plasma evolution of the 
foam sample is studied and the implications for the ion- 
stopping experiment are discussed. 

2. RALEF-2D 

2.1. Two-dimensional radiation hydrodynamics 

All presented results have been obtained with a newly 
developed radiation-hydrodynamics code RALEF-2D 
(Radiative Arbitrary Lagrangian-Eulerian Fluid dynam- 
ics in two Dimensions) 11811 . whose hydrodynamics part 
is based on an updated version of the CAVEAT hydro- 
dynamics package 11911 . The one-fluid one-temperature 
hydrodynamic equations are solved in two spatial di- 
mensions (in either Cartesian (x,y) or axisymmetric 
(r,z) coordinates) on a multi-block structured quadri- 
lateral grid by a second-order Godunov-type numerical 



scheme. Mesh rezoning and remapping is performed 
within the Arbitrary Lagrangian-Eulerian (ALE) ap- 
proach to numerical hydrodynamics. 

Thermal conduction and radiation transport have 
been implemented by newly developed algorithms 
within the unified symmetric semi-implicit approach 
112011 with respect to time discretization. For the ther- 
mal conduction, a conservative, second-order accurate 
symmetric scheme on a 9-point stencil 12 ill is used. For 
the radiation energy transport, the quasi-static transfer 
equation 



Q..VI v = k v (B v -I v ) 



(1) 



for the spectral radiation intensity I v - I v (t,x,£l) is 
solved numerically in order to couple the radiative heat- 
ing term 
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to the hydrodynamic energy equation. Within the quasi- 
static approximation the term c~ l dl v /dt (where c is the 
speed of light) is omitted in equation (Q]). Spatial dis- 
cretization and integration of equation (JTJ is achieved by 
the classical S„ method 12211 to treat the angular depen- 
dence of the radiation intensity I v (t, x,il), and by the 



method of short characteristics I123ll . which ensures that 
every grid cell automatically receives the same number 
of light rays. The correct reproduction of the diffusion 
limit on distorted non-orthogonal grids I124ll is guaran- 
teed by a special combination of the first- and second- 
order interpolation schemes in the finite-difference ap- 
proximations to equations ([TJ and (0. 

Energy deposition by a monochromatic laser beam is 
described by means of the inverse bremsstrahlung ab- 
sorption. Numerically, propagation of the laser light is 
treated within the same algorithm as the radiation en- 
ergy transport — which means without refraction. 

2.2. Equation of state and opacities 

The equation of state, thermal conductivity, and spec- 
tral opacities used in the present work were provided 
by the THERMOS code lH, which has been devel- 
oped at the Keldysh Institute of Applied Mathematics 
(Moscow). The spectral opacities are generated by solv- 
ing the Hartree-Fock-Slater equations for plasma ions 
under the assumption of equilibrium level population. 
In combination with the equilibrium Planckian intensity 
By, used in (Q]) as the source function, the latter means 
that the radiation transport is treated in the local thermo- 
dynamic equilibrium (LTE) approximation. The appli- 
cability of the LTE approach to our problem is discussed 
below in section [4] 
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Figure 1: Spectral absorption coefficient k v of gold at T = 100 eV and 
p = 10 mg cm used in the simulations: the original THERMOS 
code data (thin solid curve) are shown together with the group- 
averaged values for 24 (thick solid curve) selected spectral groups. 



The transfer equation ([TJ is solved numerically for 
a selected number of discrete spectral groups \vj, Vj+i , 
with the original THERMOS absorption coefficients k v 
averaged inside each group j by using the Planckian 
weight function. Two different sets of frequency groups 
are prepared for each code run: the primary set with 
a smaller number of groups (24 in the present simula- 
tions) is used at every time step in a joint loop with 
the hydrodynamic module, while the secondary (diag- 
nostics) set with a larger number of groups (300 in 
the present simulations) is used in the post-processor 
regime at selected times to generate the spectral out- 
put data. Figure [TJ gives an example of the spectral de- 
pendence of k v for a gold plasma at T - 100 eV and 
p = 0.01 g cirT 3 , together with the 24 group-averaged 
values. 



3. Target geometry and simulation setup 

3.1. Three-dimensional experimental target geometry 

The three-dimensional (3D) configuration of the 
combined hohlraum-foam target used in the experi- 
ments [17!] is shown in Fig. [2] The hohlraum wall 
is a thin cylindrical gold layer with the inner diame- 
ter dh r = 1.3 mm and the length Z/„. = 2.0 mm, em- 
bedded inside a massive aluminum block. A 1.4-ns 
long frequency-doubled (Ai — 526.5 //m) PHELIX laser 
pulse with 0.2 ns long flanks (Figure[3]l shoots at an an- 
gle of a/ = 45° with the total pulse energy of E\ - 180 J 
to the center of the hohlraum wall. The spatial profile 
of the laser intensity was approximated by a Gaussian 
curve with a full width at half maximum of 0.2 mm. 




Figure 2: (color online) 3D scheme of the hohlraum-foam target: 
shown is a cylindrical gold hohlraum embedded into an aluminum 
block; the PHELIX laser (green) shoots from the top; a cylindrical 
cellulose-triacetate (C^HigOg) foam at the bottom of the hohlraum 
is surrounded by a copper holder; after creation of the hot plasma, the 
stopping power for the ion beam (blue dashes) is measured. 



At the bottom of the hohlraum, a cylindrical sam- 
ple made of a cellulose-triacetate (CnH^Os) foam with 
initial mean density (pj) - 2.0 mg/cm 3 , diameter 
df — 1.5 mm, and length If - dj lr = 1.3 mm is attached. 
Except for the entrance and exit holes for the ion beam 
(blue dashes in the figure) and the connecting hole to 
the hohlraum, the foam is surrounded by a solid copper 
holder. The chosen foam density is, on the one hand, 
sufficiently low to ensure efficient propagation of the 
hohlraum radiation, and, on the other hand, high enough 
to allow fabrication of stable 3D samples by tempera- 
tures up to 220 °C; the finest pore structure within the 
foams cellular network is on the order of 1 yum. Note, 
that the column density of the used foam configuration 
corresponds to that of a 1.3 /mi thick solid carbon foil, 
but the combination of the low density and the large 
thickness of the foam provides the advantage of a di- 
minished hydrodynamic expansion on the nanosecond 
time scale in comparison to the foil. 

In the experiments, the ion bunches of duration 
tf, — 3 ns and diameter db - 0.5 mm probed the hot 
plasma, and the heavy-ion energy losses in the ionized 
sample measured by the time-of-flight method were 
compared with the previous ion shot through the cold 
foam. Simple evaluation of the corresponding Coulomb 
logarithms indicates that, for the 4.77-MeV/u titanium 
ions used in the experiments, the stopping power of the 
hot foam (approximated as 3.75-times ionized pure car- 
bon at p - 2.0 mg/cm 3 and T ~ 25 eV) should exceed 
the corresponding cold value by some 75%. 
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Figure 3: Temporal profile of the deposited laser power W c i ep for the 
simulated two-dimensional configuration. 



3.2. Two-dimensional simulation setup 

Theoretical modeling of the above described target 
with an oblique (with respect to the hohlraum axis) in- 
cidence of the laser beam is intrinsically a 3D problem. 
Having no 3D radiation-hydrodynamics code at hand, 
we had to reduce it to two dimensions. Figure |4] shows 
the corresponding lateral cut of the actually simulated 
2D target. 

Numerical simulations were performed in the {x,y) 
Cartesian coordinates, where the 3D extension of the 
simulated region spreads to infinity along the z-axis. 
The 2D hohlraum wall is represented by two gold plates, 
each with the initial density pa u - 18 g/cm 3 . The 
cellulose-triacetate foam was modeled as a homoge- 
neous block of pure carbon with the initial foam density 
of pc — 2 mg/cm 3 , supported from below by a copper 
plate with pc u = 8 g/cm 3 . Two additional horizontal 
gold plates are added to confine the lateral plasma ex- 
pansion. All the dimensions indicated in Fig.|4]coincide 
with the 3D experimental values. The supposedly empty 
parts of the simulated domain were initially filled with 
corresponding gases at sufficiently low densities, so that 
their overall dynamic and thermal influence were negli- 
gible. Simulations started from the initial state of pres- 
sure equilibrium among all the target parts; the bound- 
ary condition of free outflow was applied at all the outer 
edges. 

Spectral energy transport by thermal radiation was 
treated with 24 logarithmically equidistant (except for 
the first group) discrete frequency groups delimited by 
the photon energies 



hvj = 1(T 4 , 0.02, 0.026, ..., 10.0 keV, 



(3) 



as is shown in Fig.Q] The delimiting frequencies for the 
300 diagnostics groups were logarithmically uniformly 
distributed between 0.01 keV and 10 keV. The angu- 
lar dependence of the radiation intensity was calculated 
with the S 30 method, which offers 960 discrete ray di- 
rections over the entire 4n solid angle. 



/ 




AufillCIO^gcrri ) 
Au (18 gem" 3 ) 
(x,y) = (0,0) 




cm" ) 



cm" ) 



Cu (8 g cm"* 



Figure 4: (color online) Lateral cut of the simulated 2D configuration; 
all dimensions and initial densities are given; the whole configuration 
extends to infinity along the z-axis; the hohlraum is represented by two 
gold plates, the foam by a carbon block supported by a copper plate; 
for simplicity two horizontal gold (rather than aluminum) plates are 
used to ensure the upper confinement of the lateral plasma expansion. 



3.3. Reseating of the laser input energy from 3D to 2D 

In 2D simulations we used the same spatial and tem- 
poral laser power profiles as in the experiments, and 
only the total input energy was rescaled to ensure that 
in two dimensions the hohlraum walls absorb approxi- 
mately the same amount of energy per unit surface area 
and exhibit similar dynamics of wall evaporation as in 
the original 3D configuration. For the 3D input energy 
of Ei - 180 J our rescaling procedure, described in 
|Appendix A| yields (after one iteration) a 2D input en- 
ergy of EY' = 122.8 J/mm, which was used in all 2D 
simulations. The 2D runs needed to iterate Ei were per- 
formed for the isolated hohlraum without the attached 
carbon block. The corresponding areas of the hohlraum 
walls and holes, entering the equations (!A.4t and ( 1A.7I >. 
are 



S w - ndhrhr, 
o w — llhr, 



Sh - 2dh r . 



(4) 
(5) 



In Section [4] we present the results for the isolated 
hohlraum obtained with the first-order input energy E) 
on a rectangular mesh with ^ 265000 cells in the purely 
Eulerian mode. The simulation which included the car- 
bon plasma was performed in the ALE mode with La- 
grangian interfaces at the material boundaries on a mesh 
of ^ 1 10000 cells and is discussed in Section[5] 
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Figure 5: (color online) Color contour plots of a) the deposited laser power q (semi-transparent plot superimposed on the density), b) the matter 
temperature T, and c) the radiation temperature T r inside the hohlraum by the end of the laser pulse at t = 1.2 ns; by this time the expanding hot 
plasma plume at the laser spot on the left wall is heated up to the maximum temperatures T max = 324 eV and T rmax = 114 eV. 
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Figure 6: Calculated x-ray spectrum as would have been observed 
through the lower hohlraum hole and as seen by the foam sample at 
t = 1.2 and 6 ns; the dashed line shows a Planckian fit for T = 28 eV. 
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Figure 7: Temporal evolution of the average matter temperature T and 
radiation temperature T, at the center of the hohlraum (within a box 
of size 0.2 X 0.2 mm) on a time scale of 18 ns. 



4. Radiative properties of the isolated hohlraum 

The heating of the carbon plasma is determined by 
the thermal history of the hohlraum, especially during 
the laser pulse. Figure |5ji demonstrates how the laser 
light is absorbed near the end of the laser pulse, at 
t = 1.2 ns; Figs. |5p and [5}; show the matter and radi- 
ation temperatures inside the hohlraum at this moment. 
A plume of hot laser-ablated gold plasma rises from the 
left wall with an average velocity of w 3.5 x 10 7 cm/s. 
Also, some plasma is evaporated from the upper cor- 
ner of the right wall licked by the periphery of the laser 
beam. So long as this plasma remains undercritical and 
does not significantly block the laser light, its influence 
can be neglected. 



The temperature inside the expanding plasma cloud 
reaches its peak value T » 0.38 keV at t — 0.2 ns, 
the density varies in the range p*2x 10~ 4 - 5 x 10~ 3 
g/cm 3 , and the LTE ionization degree amounts to 
Z as 45 - 55, which implies a free electron density in 
the range of n e « 10 19 - 10 21 cm 4 . Recall that the 
critical electron density for a A — 526.5 nm laser light is 
n e ,cr ~ 4 x 10 21 cirT 3 . ForZ = 40, the critical density of 
free electrons in a gold plasma corresponds to the mass 
density of p cr — 3.3 x 10~ 2 g/cm 3 . 

The calculated x-ray spectra, as would have been ob- 
served through the lower hohlraum hole at t - 1.2 
and 6 ns, are shown together with a Planckian fit of 
the second spectrum for T — 28 eV in Fig. [6] This 
figure displays the spectral power F v [TW cirT 1 sr~ : 
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Figure 8: (color online) a) Color contour plot of the matter tempera- 
ture T at t = 1.2 ns, and b) vertical temperature profiles through the 
carbon block center at three selected times, ( = 1.2, 2.5, and 6 ns; the 
shaded (cyan) vertical strip of width db marks the ion beam aperture. 



keV" 1 ] per unit cylinder length, obtained by integrat- 
ing the spectral intensity I v along an imaginary obser- 
vation slit (perpendicular to the y-axis) over the range 
-0.65 mm < x < +0.65 mm. Taking its origin from the 
groups of strong emission lines of gold (cf. Fig. [1}, the 
highly non-Planckian spectrum at t - 1.2 ns exhibits 
three peaks between 0.18 - 0.3 keV, 0.78 - 0.95 keV, 
and 2.5 - 2.7 keV. As will be shown in the next sec- 
tion, such a spectrum is quite favorable for the efficient 
quasi-volumetric heating of the carbon foam. 

On a longer time scale, displayed in Fig. [7] the mat- 
ter and radiation temperatures inside the hohlraum come 
close to equilibrium shortly after the laser is off, and for 
t > 3 ns stabilize at T a T r » 30 eV: at this stage the 
hohlraum gradually cools down due to the radiative en- 
ergy loss and the calculated x-ray spectrum at t = 6 ns 
approaches the Planckian shape. The average matter 
and radiation temperatures near the hohlraum center, 
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Figure 9: Spectral absorption coefficient k v of carbon at T = 10 eV 
andp = 2 mg cm used in the simulations: the original THERMOS 
code data (thin solid curve) are shown together with the group- 
averaged values for 24 (thick solid curve) selected spectral groups. 



plotted in Fig. [7] for times t < 18 ns, were calculated 
by averaging over the grid-cell masses in the central re- 
gion (-0.1 mm < x < 0.1 mm, 0.9 mm < y < 1.1 mm). 

These results allow us to assess the applicability of 
the LTE model to our problem. In general, non-LTE ef- 
fects may only become significant when (and where) the 
local radiation spectrum strongly deviates from that of a 
black-body with the local matter temperature T. In our 
case this happens only during the short period t < 1 .4 ns 
of laser illumination and only in the vicinity of the laser 
spot. As a consequence, the non-LTE physics can be ex- 
pected to affect the ionization equilibrium of gold in the 
laser focal spot and to modify the spectral details (es- 
pecially the line features at hv > 2 keV) of the primary 
hard x-ray flash, generated at t < 1 .4 ns. However, once 
we focus our attention on later times t > 2 ns when the 
radiation field everywhere (including the carbon foam 
discussed below) comes close to local equilibrium with 
matter (see Figs. |6]and[7ji, the use of the LTE approxi- 
mation is fully justified. 

At t S3 7 ns, the central region of the hohlraum is 
slightly reheated to T ~ T r x 33 eV due to the collision 
of the expanding clouds of the ablated material from the 
two hohlraum walls. Such a collision leads to the forma- 
tion of a strongly radiating shock front ||26ll with a prac- 
tically full conversion of the kinetic energy into thermal 
radiation, accompanied by a strong plasma compres- 
sion. As a result, a thin and dense filament of shock- 
compressed gold plasma — clearly visible in Figs. [TUh 
and [10b below — is formed. In our case this filament 
stays close to the hohlraum center, which agrees well 
with the experimental x-ray pinhole images showing a 
bright spot at the hohlraums center 11171 Fig. 4b]. 
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Figure 10: (color online) Color contour plots of a) the matter density p,b) the matter temperature T, andc) the LTE ionization degree Z of the whole 
experimental configuration at t = 6 and 14 ns. In the middle of the hohlraum one observes a filamentary structure, evolving from the collision of 
the ablated gold plasma flows. Material interfaces can be identified as sharp discontinuities on the Z plot. Inside the carbon block one can discern 
an almost planar shock front, propagating upwards from the copper support plate and reaching y = -1 mm at t = 13.5 ns. 



5. Dynamics of the carbon plasma 

5.7. Heating dynamics 

As may be expected, most of the energetic x-rays, 
originating from the focal laser spot plasma, enter and 
heat the carbon block already during the laser pulse. 
The distribution of the matter temperature by the end 
of the laser pulse at t — 1.2 ns is shown in Fig. [8^. At 
this time a strong spatial variation of the carbon plasma 
temperature, ranging from values below 10 eV to those 
above 55 eV, can be observed. Furthermore, some x- 
rays, mainly those which originate from the laser spot 
at the left hohlraum wall, pass through the carbon block 
at its upper corners and "shine" into the low-density car- 
bon vapor fill. 

The vertical temperature gradient across the plasma 
and its relaxation in time are illustrated by three selected 
temperature profiles through the center of the carbon 
block in Fig. [8j>. For a proper qualitative understand- 
ing of the carbon heating process one has to consider 
the spectral absorption coefficient k v of carbon, shown 
in Fig.|9]for T — 10 eV andp = 2 mg cm 4 . By compar- 
ing the plots in Figs.[6]and[9]one sees that the first spec- 
tral peak of strong emission from the hot (T =* 300 eV) 
gold plasma in the laser spot practically coincides with 
the transparency window hv - 0.1-0.3 keV of carbon 
just below its K-edge, where k v x 1 mm -1 ; the pho- 
tons with hv > 0.8 keV from the second spectral emis- 
sion peak also have k v < 1 mm -1 . In other words, for 
a large portion of the hohlraum radiation emitted dur- 
ing the laser pulse the carbon foam has an optical thick- 
ness of w 1 . The latter means that the carbon foam is 
practically instantaneously (i.e. synchronously with the 
temporal laser power profile) heated by a flash of x-rays 
from the laser spot over the entire foam volume to an 
average temperature of T ~ 30 eV, varying by about a 
factor 4 across a distance of 1 mm. At the same time, 
about 20% of the hohlraum x-ray emission, generated 
during the laser pulse, passes through the foam and di- 
rectly hits the copper support plate. 

5.2. Temperature equilibration phase 

After the laser is turned off, the hohlraum continues 
to glow in soft x-rays on a time scale of tens of nanosec- 
onds, letting out the energy accumulated in its wall ma- 
terial during the laser pulse. Figure HOl shows the 2D 
distributions of the density p, the matter temperature T, 
and the LTE ionization degree Z for the whole exper- 
imental configuration at t - 6 and 14 ns. From these 
plots some important observations can be made for the 
ion-stopping measurements. 
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Figure 11: Normalized column density cr(t,y) along the ion trajec- 
tories parallel to the x-axis (see Figs. |2]and|4} as a function of the 
y-coordinate for four selected times t = 6, 10, 14, and 18 ns. The 
positions of the copper and gold plasma fronts are easily identified 
as almost vertical walls on the left and on the right. The shock front 
(smeared due to 2D effects) in the foam, pushed by the expanding 
copper, reaches the ion beam aperture di, at t = 13.5 ns. 



Firstly, one can clearly see that the expanding carbon 
plasma pushes back the gold plasma from the hohlraum, 
which means that one should not fear the ablated gold to 
get into the way of the ion beam. Secondly, as a signifi- 
cant energy portion of the main x-ray flash, generated by 
the laser pulse, penetrates through the translucent car- 
bon foam and is absorbed by the copper holder, the re- 
sulting pressure discontinuity between copper and foam 
launches a shock wave into the foam. This shock front 
propagates towards the hohlraum and enters the ion- 
beam aperture -1.0 mm < y < -0.5 mm at t « 13.5 ns, 
i.e. late enough to perform the measurements at earlier 
times. Finally, the simulation shows that a reasonably 
homogeneous plasma volume with Z w 3.7-3.8 and 
T ~ 23-27 eV is created along the ion beam path, quite 
suitable for the ion-stopping measurements within a cer- 
tain limited time window. 

The appropriate time window [t\,t2] for the ion- 
stopping measurements can be evaluated by inspecting 
the data presented in Figs. [TTI and [121 Figure QT] shows 
the normalized column density 



aif,y) = 



I p(t,x,y)dx I p(0, x,y)dx 



(6) 



along the supposed ion trajectories parallel to the x-axis 
(see Figs. [2] and [4]i as a function of the y-coordinate 
at four selected times t — 6, 10, 14, and 18 ns. For 
high-quality ion stopping measurements, it is important 
to have the values of cr(t , y) as close to 1 as possible for 
the entire ion pulse duration t\, » 3 ns over the entire 
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Figure 12: Temporal evolution of a) the matter temperature T and b) 
the LTE ionization degree Z at the top (y = -0.5 mm) and the bottom 
(y = - 1 .0 mm) of the ion beam cross-section in the center of the foam 
block (x = 0). Figure c) shows the maximum spatial variation of the 
normalized column densities <x(f, y) within the ion beam cross-section 
and the difference between the top and bottom ionization degree val- 
ues weighted by their mean value Z = j(Z top +Z/, ottom ) as functions of 
time. A propitious time window for the ion-stopping measurements is 
marked as a shaded (yellow) vertical strip. 



beam aperture -1.0 mm < y < -0.5 mm. The eventual 
significant departures of cr(t,y) from 1 are caused by 
three effects: (i) the 2D lateral expansion of the carbon 
plasma beyond the confining walls of the copper holder, 
(ii) compression of the foam material by the shock front 
propagating from the bottom of the carbon block, and 
(iii) swelling of the gold walls near the hohlraum edges. 
Our results indicate that for t < 8 ns the maximum de- 
viations of cr(t, y) from unity within the diameter of the 
ion beam d\, do not exceed 10%. The t — 14 and 18 ns 
profiles in Fig.QT|are already significantly perturbed by 
the effects (i) and (ii). 

Figures [12] a) and b) show the temporal evolution 
of the matter temperature T and the ionization degree 
Z at two points (x,y) - (0,-0.5) and (0,-1.0) that 
are of special interest from the point of view of the 
ion-stopping measurements. Figure Q~2] c) shows the 
difference between both ionization degree values and 
the maximum spatial variation of cr(t, y) within the ion 



beam aperture. One sees that the temperature and ion- 
ization become quite uniform for t > 10 ns — which 
would be a good time for the ion bunch to arrive if not 
already significant perturbations to the <r{t , y) values oc- 
cured. 

Finally, as a result of complex interplay between the 
above discussed physical effects, the best compromise 
for the ion measurements in the analyzed configuration 
should roughly fall into the range 3 ns < t < 8 ns — 
as is indicated in Fig. [12] with a shaded vertical strip. 
Note that the upper boundary of this window can easily 
be expanded by increasing the size of the foam block: 
every extra 100 ^m of the foam size would add roughly 
1 ns to the favorable time span. 

6. Conclusion 

The results of a 2D radiation-hydrodynamic study 
of a combined hohlraum-foam target used to perform 
ion-stopping measurements in a homogeneously ion- 
ized plasma are presented. Our principal goal was to 
investigate how uniform and for how long a sample 
plasma column can be created in this type of laser- 
driven targets. We have analyzed the principal physi- 
cal processes, which determine the spatial structure and 
dynamics of the sample plasma, and demonstrated how 
their combined effect defines an optimum — to ensure 
clean ion-stopping measurements in a plasma — time 
window for accepting a =* 3-ns long bunch of fast ions. 
For the target parameters used in this work, this time 
interval turns out to be between 3 ns and 8 ns after the 
onset of the laser pulse. We find that within this win- 
dow the time and space variations of such key param- 
eters as the column mass density along the ion trajec- 
tories and the plasma ionization degree do not exceed 
±7%. On the basis of our analysis one can readily iden- 
tify the guidelines for target modifications that could 
improve the quality of the ion-stopping measurements 
in this type of laser-plasma targets. 
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Appendix A. Reduction of a 3D hohlraum to a 2D 
configuration 

Reduction of an intrinsically 3D problem to two di- 
mensions requires an additional spatial symmetry to be 



imposed on the original 3D configuration, which may 
be either a translational invariance (along the z-axis) or 
an invariance with respect to rotation around a fixed z- 
axis. In result one would want the new 2D configura- 
tion to reproduce as closely as possible all the main dy- 
namic features of the original 3D problem. In our case 
it means that the simulated 2D configuration must rep- 
resent a certain 2D cut of the original 3D hohlraum with 
all the key dimensions taken from the latter. Having 
preserved the temporal and spatial shapes of the driving 
laser pulse, we are left with only one undetermined pa- 
rameter for the 2D case, namely, the total input laser en- 
ergy, which we will denote as Ei for the 2D case, and as 
Ei for the original 3D case. Note that the original value 
Ei must be recalculated (rescaled) to the "equivalent" 
2D value E\ already because in the case of translational 
invariance these two quantities have different physical 
dimensions. 

Since in hohlraum-type targets both the matter mo- 
tion and the intensity (as well as the spectrum) of 
the hohlraum radiation are controlled by the inci- 
dent radiation-energy fluxes per unit surface area of 
the hohlraum interior, we assume that the "physically 
equivalent" 2D hohlraum must on average accept the 
same amount of energy per unit inner surface area as 
the original 3D hohlraum. Mathematically this condi- 
tion can be expressed as the following two equations of 
the global energy balance in a hohlraum 



Ei - F W S W + FhSh, 
Ei = F W S W + FhSh- 



(A.l) 
(A.2) 



Since the surface areas S w , S/,, S w , and Si, are all 
known, Eqs. (1A.U and ( IA.21 > yield the following ex- 
pression for the rescaled input energy 



E, = E 



Sw + qimSh 

S w + qim-Si, 



where 



Clin 



Fh 

F w 



(A.4) 



(A.5) 



is an unknown dimensionless factor. It can be found by 
the method of successive approximations as follows. In 
the zeroth approximation, one can set F/, - F w , i.e. as- 



sume that qi n 
with 



<lh 



(0) 



1 , and perform a 2D simulation 



E, = m = Ei S -^±. 
1 s H , + s h 



(A.6) 



Making use of the results of this simulation, one can 
evaluate the first-order value of q/ m , as 
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(i) 



e(0) 

*" r.out 



§ h £f> 



#(0) 



(A.7) 



e(0) 



where E rout is the total amount of radiative energy, 
which escapes the hohlraum in the zero-order 2D run. 
Equation ( IA.7t is a direct consequence of the global en- 
ergy balance relation (I A. 2\ combined with 



t^r,out — r h& I 



(A.8) 



Here S w and Sh are, respectively, the total surface ar- 
eas of the inner 3D hohlraum wall and of all its holes, 
S w and § i, are the corresponding quantities in the 2D 
hohlraum, F w is the radiation-energy fluence (measured 
in J/cm 2 ) absorbed by the inner hohlraum walls, and 
F}, is the radiation-energy fluence which escapes the 
hohlraum through its holes. If the hohlraum radiation 
is characterized by the equivalent black-body tempera- 
ture T r (t), then F/, can be calculated as 



A. 41). one obtains a 



Fi, = o~sb 



J T ; 



ell. 



(A3) 



where osb is the Stefan-Boltzmann constant. As a word 
of caution it should be noted that, because it is only 
the time- and space-averaged quantities F w and F/, that 
are ensured to have the same values in the original 3D 
hohlraum and its 2D counterpart, one cannot expect that 
all the details of the 3D problem should be adequately 
reproduced by a corresponding 2D simulation. 



Having substituted q;' into Eq. 
first-order estimate E { t ' for the rescaled input energy £/. 
Higher-order approximations are usually not needed. 
Note that hohlraums with walls of a heavy metal typ- 
ically have qi m , w 2-3. 
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